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I. DEFINITION OF THE SUBJECT AND ITS IMPORTANCE 

The terminology granular matter refers to systems with a large number of hard objects 
(grains) of mesoscopic size ranging from millimeters to meters. Geological examples include 
desert sand and the rocks of a landslide. But the scope of such systems is much broader, 
including powders and snow, edible products such a seeds and salt, medical products like 
pills, and extraterrestrial systems such as the surface regolith of Mars and the rings of 
Saturn. The importance of a fundamental understanding for granular matter properties can 
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hardly be overestimated. Practical issues of current concern range from disaster mitigation 
of avalanches and explosions of grain silos to immense economic consequences within the 
pharmaceutical industry. In addition, they are of academic and conceptual importance as 
well as examples of systems far from equilibrium. ^ 

Under many conditions of interest, granular matter flows like a normal fluid [1]. In the 
latter case such flows are accurately described by the equations of hydrodynamics. Attention 
is focused here on the possibility for a corresponding hydrodynamic description of granular 
flows. The tools of nonequilibrium statistical mechanics developed over the past fifty 



years for fluids composed of atoms and molecules 0, S], are applied here to a system of grains 
for a fundamental approach to both qualitative questions and practical quantitative predic- 
tions. Applications of basic atomic physics principles to granular fluids have accelerated 
during the past decade, starting with an emphasis on molecular dynamics (MD) simulations 



61 and kinetic theory 

nnl 

described here [9|, HO, 111 



, and more recently with the theoretical methods of the type 
ld,[l3|,Q. 



II. INTRODUCTION 



To start with the familiar, consider a jar of vitamin pills, mustard seeds, or peanuts. 
Remove the lid and pour them into a bowl, observing that the "flow", or their collective 
motion, has some similarity to that of a normal fluid such as water. The collective motion 
in both cases is the consequence of collisions among their constituents, grains or atoms, 
and their large number. It is tempting to make the correspondence of grains to atoms 
in considering the similarities of flows in these two types of fluids. The objective here is 
to explore in formal detail the extent to which that correspondence is conceptually and 
quantitatively justified. An important prerequisite is the integrity of the grains during their 
motion. Each grain is comprised of a large number of atoms or molecules. Integrity refers 
to their retention of mass and shape following interactions with other grains or with their 
environment. As such, the grains behave as "particles" whose detailed internal structure is 
not essential to their description, which is captured instead by a few parameters describing 
their shape, mass, and collisional properties with other grains. However, an important 
consequence of this underlying molecular structure is a redistribution of translational kinetic 
energy of the grains and internal energy of the constituent molecules. At the mesoscopic 
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level this appears as an energy loss on collisions between pairs of grains. This is a central 
feature of granular fluids differentiating them from atomic fluids: the inelasticity of granular 
pair collisions. 

Granular matter occurs in two classes of states, compact and activated In the 

first case, the grains form a static packed configuration within the container due to the 
effects of gravity on their relatively large mass and their inelastic collisions. Any initial 
motion is quickly dissipated and their kinetic energy becomes negligible relative to the 
gravitational potential energy. Important questions arise about the possible and probable 
packing configurations that determine the stresses within the system and the distribution 
of forces on the container. For example, chains of particles in contact can occur as arches 
to support matter above them while reducing their force on the matter below. There is an 
intense interest in the study of such states, generically referred to as contact mechanics. 

Activated states refer to continuously driven systems, or gravity free conditions. For 
example, a container of grains in a compact configuration can be shaken to impose kinetic 
energy and motion among the grains. Similarly, unrestrained systems in a gravitational 
field will flow towards lower potential energy (e.g., hopper flow or flow down an incline). 
Initial activation in space laboratory experiments provides another example (self-sustained 
fluidization). For the flows considered here as candidates for a hydrodynamic description 
continual collisions are essential. This means that within each small cell, still containing 
many particles, the particles are moving randomly relative to the collective motion of that 
cell. Thus, ballistic motion or beams with all particles moving independently in the same 
direction are excluded. 

Both compact and activated grains may occur immersed in a continuum such as water or 
air that may have a strong or weak effect on their collective properties. For compact systems 
water may provide a lubrication effect that affects the dominant class of configurations. For 
activated systems it can provide an additional dissipative drag between collisions among 
the grains. When the medium plays an important role the systems is said to be wet. In 
the opposite limit it is said to be dry. Finally, it is possible for compact and activated 
components of a system to coexist as heterogeneous states. Here, only the simplest case of 
dry systems in fully activated flows are considered. These are referred to in the following as 
granular fluids. 

Advances in the study of granular fluids have arisen from many communities, including 
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chemical engineering, materials sciences, and physics. The additional academic and con- 
ceptual importance of granular matter as practical systems for exploring the relevance of 
many-body fluid methods is primarily for the physics community. Granular matter, viewed 
as a system of particles with inelastic interactions, provides new opportunities to test the 
qualitative and quantitative limits of many-body methods developed over the past century 
: or atomic and molecular systems. This is the field of non-equilibrium statistical mechanics 
s], y, [5]. Granular matter provides a new testing ground for a reconsideration of the most 
fundamental concepts and tools with the potential for enhanced understanding of their 
place in atomic and molecular systems as well. 

Statistical mechanics addresses the difficult many-body problem of extracting macro- 
scopic properties of experimental interest from the very large number of constituent parti- 
cles. The results express these properties in terms of the fundamental "microscopic" features 
of these particles, such as mass, shape, degree of inelasticity, and collisional properties. In 
this way the properties of the vitamin pills, mustard seeds, and peanuts are distinguished 
at a fundamental level. Also, conceptual issues such as the limitations of a macroscopic de- 
scription are exposed through the insistence on their logical evolution from the fundamental 
microdynamics. In the next section, the granular fluid is described as a system of particles 
interacting via pairwise additive, nonconservative forces. The microscopic dynamics of these 
particles leads to balance equations for the mass density, energy density, and momentum 
density. Their averages define the "hydrodynamic fields" which are candidates for a macro- 
scopic, continuum mechanics description. These exact equations are described in Section 
3 and the need for "constitutive equations" to provide a closure is described. The origin 
of constitutive equations, and consequently the origin of hydrodynamics, is associated with 
the concept of "normal states" in Section 4. The normal state for the case of small 
spatial deviations from homogeneity is constructed formally in Section 5, resulting in the 
constitutive equations for Navier-Stokes hydrodynamics [17\. This derivation also provides 
insight into the context in which such a description should hold, and differences from the 



corresponding results for a normal fluid are noted. Empirical evidence 
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22| , and corresponding results from kinetic theory 



3, 
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25 



19| . simulations 
26| support the appli- 



cability of this hydrodynamic description under appropriate conditions. Finally, the results 
are summarized in Section 6 and some comments on the outlook for future developments 
are offered. 
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The presentation here is focused on recent work of the author and his collaborators for 
application of statistical mechanics to explore hydrodynamics for a granular gas. Conse- 
quently, the references quoted are heavily weighted toward those developmental studies. 
Apologies are offered for the exclusion of the vast and important complementary literature 
on simulations, kinetic theory, and experiments also bearing on this topic. Many of these 
can be found in the list of Books and Reviews given here. 

III. GRANULAR FLUID AND ITS STATISTICAL MECHANICS 
A. Nonequilibrium statistical mechanics 

Consider a system of >> 1 identical grains (hereafter referred to as particles) in a 
volume V, whose initial positions {qj} and velocities {vj}, 1 < i < N, are specified. The 
positions and velocities define a point in a 6A^ dimensional space denoted by F = {qi,Vj}, 
defining the microstate of the system. A macrostate is defined by a probability density 
p (r) in this space, representing statistical rather than precise knowledge of the system. The 
field of statistical mechanics addresses properties of macrostates, based on the recognition 
that for very large the details of microstates are neither experimentally accessible nor 
practically calculable. Properties of interest are represented by functions ^(r), and their 
values for a macrostate p (F) are determined from the expectations 

{A;p)^ J dTpiT)AiT). (3.1) 

In this section, a brief overview of the essential ingredients of nonequilibrium statistical 
mechanics is given, broadened from its usual form [3| to include granular matter. 

The dynamics of macrostates is determined from the underlying dynamics of the mi- 
crostates. The initial point F changes in time since the particles have velocities and move to 
new positions. They move in straight lines until one or more come within the force field of 
other particles, at which point their velocities change as well as their positions. The forces 
are taken to be pairwise additive, such that the total force on particle i is Fj = Fij, where 
Fij is the force on particle i due to particle j. This does not mean that the interactions are 
pairwise sequential; three or more particles can interact simultaneously. The pair forces are 
restricted by Newton's third law, Fij = —Fji, with conservation of momentum. Otherwise 
quite general forces can be considered to represent the shape of the particles and their degree 
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of inelasticity. It is assumed here that the force range vanishes outside a distance (t/2 from 
the center of each particle so that a characterizes the size of the particles. Furthermore, 
the particles are taken to be strongly repulsive so that their mean maximum overlap d on 
collision is small compared to their size, d/a <1. However, their size can be large or small 
compared to the mean distance between particles (V/iV)^^^, depending on whether the den- 
sity of the system is small or large, respectively. Most importantly for the purposes here, 
these forces do not conserve energy. This property captures the feature of real grains that 
center of mass kinetic energy is lost as they distort during pair collisions. Further details of 
the force law are not required at this point. 

The dynamics consists of straight line motion along the direction of the velocity at time t 
(free streaming), until the force range of any pair of particles, say z, j, overlaps. The relative 
velocity gjj = v,; — of that pair changes according to Newton's second law for the chosen 
force law Fjj. Subsequently, all particles continue to stream freely until another pair has a 
force range of overlap, and the collisional change is repeated for that pair. In this way a 
trajectory = {qi(t), . . . , qAr(t), vi(t), . . . , VAr(t)} is generated for t > 0. This traiectory is 



unique and invertible. The statistical mechanics for a fluid of inelastic particles 

Mm 

is comprised of the dynamics just described, a macrostate specified in terms of a probability 
density p(F), and a set of observables generically denoted by A{T). The expectation value 
for an observable at time t > for a state p(F) given at t = is defined by 

(A(t);0) = j dVp{V)A{Vt) = j f/Fp(F)e*M(F) (3.2) 

where A{t) = A(Tt), and Ff = {qi(t), . . . , qjv(t), vi(t), . . . , vjv(t)} is the phase point evolved 
to time t from F = F(=o- The dynamics can be represented in terms of a generator L defined 
by the second equality of (13. 2p . There are two components to the generator, corresponding 
to the two steps of free streaming and velocity changes due to interactions 

N ^ N N 

L = J2^^-^^ + ^Y.Y.^^r ( Vv. - V.J . (3.3) 

i=l i=l jj^i 

An alternative equivalent representation of the dynamics is obtained by transferring the 
dynamics from the observable A{T) to the state p(F) by the definition 

J dTp{T)e'^A{T) = j (e"*V(r)) ^(H = j dV p{T ,t)A{T). (3.4) 
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The representation in terms of a dynamical state p(r, t) is referred to as Liouville dynamics. 
Its generator L is the formal adjoint of L which is found to be 

^ N N 

^ = ^ + ^ E E (^v. - V.,) • F,, (3.5) 

The difference between L and L arises because the forces are non-conservative and therefore 
depend on the relative velocities of each pair as well as there positions. Time correlation 
functions for two observables A and B are defined in a similar way 

{A{t)B; 0)^JdT (e*M(r)) p(r)i?(r) = J dTAiT) (e-*^p(r)) (e^*^i?(r)) . (3.6) 

or 

{A{t)B-0) = {AB{-ty,t) (3.7) 

In summary, averages like {A{t);0) and correlation functions {A{t)B;0) are the central 
properties of interest for a macroscopic description of physical systems. The microscopic 
dynamics can be represented in terms of the observables A{r, t) or the states p(r, t) which 
are determined from specified initial values and the equations 

{dt - L) A{V, t) = 0, {dt + L) p(r, t) = 0. (3.8) 

In the following most of the analysis is done in terms of the of states, and the associated 
equation of motion is known as the Liouville equation. 

B. Liouville equation and cooling 

For an isolated system, the total energy decreases monotonically due to the loss of energy 
on each pair collision. This is reflected in a decrease of the average kinetic energy of the 
particles between collisions and hence is referred to as coUisional " cooling" . The energy per 
particle at time t and its loss are 

e{t) = N-^{E;t), uj{t) = -dte{t) = N-^{LE;t). (3.9) 

This cooling effect is common to all solutions to the Liouville equation and it is useful to 

separate the dynamics into that due to this cooling and the residual time dependence 

pir,t) = p{r,eit),t). (3.10) 
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The Liouville equation then can be written 

dtpiV, e, t) I + {-u (e, t) + I) p(r, e, t) = 0. (3.11) 

The time derivative is now taken at constant e. The notation uj (e, t) reflects the fact that it 
is a linear functional of p(r, e, t), from its definition (13. 9p . This is a useful form that isolates 
a primary effect of the nonconservative forces (cooling) from the residual dynamics that will 
be associated with relaxation of the spatial inhomogeneities of interest below. For notational 
simplicity (13.111) is written 

((9t + Z)p(r,e,t) = 0, Z = -cj (e, t) (9e + I. (3.12) 
The corresponding equation for observables is 

(9t - £) A(r,e,t) = 0, C = -d,uj{e,t) + L. (3.13) 
where it is understood that operates on everything to its right. 

C. Stationary homogeneous state 

An isolated normal fluid supports an equilibrium state. This is a stationary solution to 
the Liouville equation with translational invariance, the Gibbs states. From the discussion 
above it is clear that isolated granular fluids have no truly stationary state due to cooling. 
However, there is a "universal" homogeneous state similar to the Gibbs state in the sense 
that a wide class of homogeneous initial states rapidly approach this state, on the time scale 
of a few collisions per particle. It is simple in the sense that all of its time dependence is 
that associated with cooling 

Po(r,t) =po({q.„v.},e(t)). (3.14) 

Here, = — q^ so the solution also has translational invariance. In the representation 
(I3.12P it is seen to be a stationary solution to the Liouville equation 

£po = 0. (3.15) 

There is no longer any explicit time dependence since uj (e, t) = oj (e) and C = —to (e) 0^ + L 
for this state. This solution is referred to as the homogeneous cooling state (HCS). Clearly, 
it is the close analogue of the Gibbs state for a normal fluid. It is an example of a "normal" 
state in the sense that all of its time dependence occurs through one of the hydrodynamic 
flelds (the energy). This concept is sharpened below. 
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IV. MACROSCOPIC BALANCE EQUATIONS 



The origins of a macroscopic description for a fluid are the balance equations for the 
average mass density (m (r) ; t), energy density (e (r) ; t), and momentum density (g (r) ; t), 
where r denotes an arbitrary field point within the system [l^ . These will be referred to as 
the hydrodynamic fields since they are the ones that are expected to obey the hydrodynamic 
equations under appropriate conditions. The phase functions m (F, r) , e (F, r) , and g (F, r) 
are well known and their explicit forms will not be needed here. They will be denoted 
collectively by Oq, (r) 

a„(r,t) <^ {m(r) ,e(r),g(r)} . (4.1) 



It follows from (13. 8p that they obey the microscopic balance equations 14 1 



dtttair, t) = Lttair, t) = -V ■ b„(r, t) - 5a2w{T, t). (4.2) 

To obtain this result, it has been recognized that the quantity Laa{r,t) can be written as 
the sum of a divergence V ■ ha{r,t) plus a remainder w{r,t) that cannot be so represented. 
For a normal fluid w{r,t) vanishes and (14. 2 p become the local conservation laws for mass, 
energy, and momentum. The bQ,(r, t) are the corresponding fluxes. This clarifies why aQ,(r, t) 
are selected for a macroscopic description. Their time dependence is determined by the 
scale of the spatial gradients of the fluxes, and averages of the latter become small as the 
system approaches homogeneity. Consequently, on long time scales the {aa (r) ; t) are the 
only surviving dynamical variables, and it is under these conditions that these fields obey 
hydrodynamic equations. The mass and momentum are conserved for a granular fluid as 
well, but there is a loss of energy w{r,t) due to the non-conservative forces. It is no longer 
obvious that the energy is still one of the slow variables since its time scale is coupled to 
w{r,t) which does not become small for nearly homogeneous states. Thus, an additional 
requirement for the existence of a macroscopic description in terms of {aa (r) ; t) is that the 
time scale of (e (r) \ t) / {w (r) ; t) must be larger than that for non-hydrodynamic properties. 
This issue is discussed further below. 

The macroscopic balance equations follow from the averages of (14.21) 

dt (a, (r) ; t) + V ■ (b, (r) ; t) = -6^2 {w (r) ; t) . (4.3) 

These equations are formally exact, but of little practical use as they do not form a closed 
(self-deterministic) set of equations for {aa{r)]t). Closure requires expressing the aver- 



age flux {ha (r) ; t) and energy loss {w (r) ; t) as functionals of the flelds {a^ (r) ; t) . Such 
relationships are called "constitutive equations". The combination of the exact balance 
equations with some form of constitutive equations provides the most general deflnition of 
hydrodynamics. 

Construction of the constitutive equations is simplifled by extracting the effects of con- 
vection. The velocity U(r, t) of a cell at point r is deflned in terms of the average momentum 

(g(r);t) = (m(r);t)U(r,t). (4.4) 

The fluxes are functions of the positions and velocities (r) = (r;{qj, v^}) = 
ha (r;{qj, Vj + U(r)}), where the velocity in the local rest frame has been introduced, Vj = 
Vj — U(r, t). Then deflning the microscopic flux in the rest frame by b^ (r) = b^ (r;{qj, Vj}) 
it follows that the average flux has the form 

(b„ (r) ; t) = (b; (r) ; t) + c,,U(r, t) ■ (b; (r) ; t) + U(r, t)da ({(a. (r) ; t)}) . (4.5) 

The flrst term is the flux of mass, energy, and momentum in a fluid element at rest, and 
represents the dissipative processes. The second and third terms are proportional to the 
flow velocity U(r,t) and are associated with convection. The coefficients of these terms 
are explicit functions of the flelds (oq, (r) ; t) (as is U(r,t)). For a normal fluid, neglect 
of the rest frame fluxes leads to the perfect fluid Euler hydrodynamic equations. Hence, 
determination of the constitutive equations is reduced to expressing the rest frame fluxes 
and energy loss as functionals of the flelds. 



V. "NORMAL" STATES AND HYDRODYNAMICS 

A hydrodynamic description is a closed set of equations for the hydrodynamic fields, 
{aa (r) ; t) . This follows from the exact macroscopic balance equations if the energy loss and 
fluxes can be represented as functionals of these flelds 

{w (r) ■,t) ^ uj{r \ {aa;t)), (b^ (r) ; t) ^ /^^(r | (a^; t)). (5.1) 

The arrow is used to indicate that such a functional representation need not be valid on 
all length and time scales, and any such restrictions constitute the domain of validity for 
hydrodynamics. The notation here and below is such that / (r, t, {(ac (r) ; t)}) denotes a 
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function of r,t and of the fields (a^ (r) ; t) at the point r, while / (r,t | {aa',t)) denotes a 
function of r,t and a functional of the {aa',t) at all space points. With such constitutive 
relations the macroscopic balance equations (14.31) become hydrodynamic equations 

dt (a« (r) ; t) + V ■ /5(r | (a«; t)) = -Sa2Uj{r \ (a,; t)). (5.2) 

The average energy loss and fluxes are averages of specific functions of the particle posi- 
tions and velocities, and hence are linear functionals of the solution to the Liouville equation. 
The existence of constitutive equations is therefore related to a special property of the so- 
lution which will be called "normal" (this terminology originates in a related context for 

n 

derivation of hydrodynamics from the Boltzmann kinetic equation [3[). The class of "nor- 
mal" distributions is defined by the functional forms 

pn (r, t) = pn ({%, Vj} I (a^; t)) . (5.3) 

All time dependence and all the breaking of translational invariance for normal states occurs 
only through the hydrodynamic fields. A familiar example of a normal distribution for real 
fluids is the local Gibbs distribution 



peiiV I (a„;t)) =exp!^q- j dr?/« (r, t) a^, (r)| 



(5.4) 



Here g is a normalization constant, and ya (r,^) are conjugate fields determined by the re- 
quirement that the averages of aa (r) give the specified values {a^ (r) ; t). In this way ya (r, t) 
are functionals of the hydrodynamic fields and pee (r | {aa', t)) is normal. The importance of 
normal solutions is that they yield directly the desired functionals of (15. ip 



u{r I [aa] ^)) = j ^^Pn (r I [aa] t)) w{r) (5.5) 

/3„(r I {aa; = j ^^Pn (r I t)) ha (r) (5.6) 

The normal state in (15.51) and (15. 6p must be a solution to the Liouville equation. In 

general, the time derivative in the Liouville equation can be separated into that which 
occurs through {aa] t) plus the residual time dependence, generalizing (I3.10p 

p(r,t)=p(r,t| {aa;t)). (5.7) 

The Liouville equation then becomes 

dtp \{a^;t) - j rfr^^-^^^{V- (b,(r);t) + 5,2(ti^(r);t)} + Ip = 0. (5.8) 
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A normal solution results when dtPn \(ac,;t)^ 0. For specified fields, (15. 8p becomes an equation 
for the r dependence of the normal phase space density as a functional of the fields. This 
dependence then allows determination of the normal forms in flS.Sp and fl5.6p . Finally, with 
the form of the hydrodynamic equation determined at that point, their solution with suitable 
initial and boundary conditions provides the explicit forms for the fields, and completes the 
normal solution. The existence and determination of this solution is the central problem for 
establishing a hydrodynamic description for both normal and granular fluids. 

The concept of a normal solution and its use in the macroscopic balance equations makes 
no special reference to whether the fluid is normal or granular, and is not restricted to 
states near homogeneity. In this general context, hydrodynamics is not a simple set of local 
partial differential equations such as the familiar Navier-Stokes equations. The latter are a 
special case of this more general idea, and their inadequacy for some conditions should not 
be interpreted as the absence of a more complex hydrodynamic description. 

In closing this Section a qualitative explanation of why a normal solution can be ex- 
pected is given, by analogy with the similar expectation for atomic fluids. For a wide class 
of initial states there is a first stage of rapid velocity relaxation in each small region toward 
the universal homogeneous state (HCS or Gibbs, respectively). However, the hydrodynamic 
parameters of that universal state are specific to each region so it is only locally homoge- 
nous, as in (15.41) for the atomic fluid. Subsequently, these differences in the parameters of 
neighboring cells are decreased by the fluxes of mass, energy, and momentum across their 
boundaries. It is this second stage where a normal description in terms of the hydrodynamic 
fields can be expected, indicating also that the space and time scales for a hydrodynamic de- 
scription should be large compared to those of the first stage. This basic conceptual picture 
is essentially the same for both atomic and granular fluids, and the rapid approach of the 
first stage is indeed observed in molecular dynamics simulation studies of both equilibrium 
and the HCS. 

VI. NAVIER-STOKES APPROXIMATION 

Equation (15. Sp presents a formidable problem and further progress requires specialization 
to specific cases of interest. Perhaps the simplest of these are weakly inhomogeneous states. 
These are states for which all spatial gradients of first order are small and all higher order 
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derivatives are negligible. Small gradients means that the relative change in the hydrody- 
namic fields over the largest microscopic length scale io is small: iodr In (a„; t) << 1. There 
are two characteristic length scales, the mean free path and the grain diameter. For a dilute 
gas the mean free path is largest, while for a dense fluid the grain size is largest. Under 
these conditions a solution to the Liouville equation can be sought as an expansion to leading 
order in these small gradients. This will be referred to as the Navier-Stokes approximation. 

According to the discussion at the end of the last section, a normal solution is expected 
after the system has relaxed to its local HCS form, denoted by po£ (r | {aa, t)), representing 
the fluid as having each cell in its own HCS. Deflne the deviations of the hydrodynamic 
fields from some common reference value by 

5 {aa;t) = {aa;t) - aoa, (6.1) 

where a^a is the same for all cells. Then, the local HCS must satisfy the conditions 

Poe (r \aoa + S {aa, t)) \s(ac,;t)=o= Po 0^; floo) , (6.2) 

dpo f , 6poi (r I ao„ + 5 (a„; t)) 

dr — 77VT7\ l<5(a„;t>=o, ■ • ■ (6.3) 



daoa J S {tta (r) ; t)„ 

i.e., the local HCS and all of its functional derivatives must agree with those of the HCS 
in the homogenous limit. Also, as a normal distribution its time dependence is through the 
exact hydrodynamic fields for the fluid state considered. This means the averages of the 
corresponding microscopic flelds (r) for the local HCS and for the solution to the Liouville 
equation must be the same 

j dT{p- Poe) a, (r) = 0. (6.4) 

A more complete discussion of the construction of poe from knowledge of po is given elsewhere 
[l^ . For the purposes here properties (16.21) . (16.31) . and (16.41) are sufficient. 

The local HCS distribution, pq£, is not a solution to the Liouville equation except in limit 
that all hydrodynamic fields become the same for each cell. Instead, it is a reference state 
approximating the actual solution after its first stage of velocity relaxation. To construct a 
solution p define its deviation from po£ by 

p(T,t\ {aa, t)) = poi (T, t I (a^; t)) + A(T,t\ {a^] t)) . (6.5) 
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The Liouville equation flS.Sp gives 

dA - [ dr'—^P-— {V ■ (b, (r) ; t) + 6^2 {w (r) ; t)} + lA 
J d {aa (r'j ; t) 

= I A^^TTwV ■ (r) ; t) + 5^2 {w (r) ; t)} - Lp,,. (6.6) 
J d (a^, (r') ; t) 

This equation is still exact, but if only small gradient states are considered it simplifies by 
retaining terms only of first order in the gradients. To be precise, the ultimate use of this 
solution is to calculate local properties of the form 

A(r, t I {y^ it)}) = J dVaiT, r)p (F, t \ (a«(r); t) +5 (a„; t)) . (6.7) 

Therefore, in the following analysis the gradient expansions are referred to the field point r 
of interest, (aQ,;t) = {aa{r);t) +6 {aa]t), i.e. the common reference values in (16. ip are the 
exact field values at the chosen point, aoa = (««(!■); t). The gradient expansion is carried out 
relative to these values. Of course the results will be general and applicable to any choice 
for r. 

The details of the gradient expansion are given in the Appendix. The solution to the 
Liouville equation to first order in the gradients is 

p{T,t\{a^;t)) = po(r,(a,(r);t)) + (l-P)(M^(r,(a,(r);t)) 

+ dt' (1 - V) (r, (a„(r); t))^ ■ V (a^ (r) ; t) . 



(6.8) 



with the definitions 



M/3 = /" dr' ( --4^1 r', (6.9) 

T„ = -(/£ + K^)^^M^. (6.10) 

The generator for the dynamics IC + has a contribution from C which is the same as in 
f l3.12p . with u evaluated for the HCS as a function of the exact hydrodynamic fields at the 
point r and time t 

C = -uJo{{aa (r) ; t))d(^e{r);t) + L. (6.11) 
The second contribution to the generator of the dynamics is the transpose of the matrix 

dujoija^ (r);t)) 
d {ap (r) ; t) 
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Finally, V is a projection operator 

dpo 



PX = ^(,J dVAfsX, Ap = j dva^ (r) , = 



d (a/3 (r) ; t) 

The phase functions Ap and form a biorthogonal set in the sense 



(6.13) 



dTA^^p = 5^p. (6.14) 

The Act are the usual global invariants of the Liouville operator L for a normal fluid; it is 
shown in the Appendix that the \E'/3 are the invariants of the new generator for dynamics in 
a granular fluid 

{ICt + K^)^P^P = Q. (6.15) 

Equation (16.81) is not quite the normal solution desired. All terms depend on time through 
{ap (r) ; t) as required, except for the last term which has an additional explicit time depen- 
dence through the upper limit of the time integral. This time dependence becomes negligible 
if the integrand is effectively non-zero after some short time scale r. Then for t » t the 
time integral becomes independent of t and can be taken formally to infinity. Thus, a normal 
solution is attained for this time scale 

Pn (r, (a„ (r) ; t)) = po (F, (a„(r); t)) + {l-V) {Mp (F, (a,(r); t)) 

+ hm r dt' (e-(^^+^")*') (1 - V) (F, (a„(r); t))] ■ V {ap (r) ; t) . 

to-^oo Jo \ / J 

(6.16) 

It is expected that the integrand should have this property of a short time scale since the 
domain of operation for the generator of time dependence is functions with translational 
invariance (as a consequence of the gradient expansion). Hence there are no explicit slow 
hydrodynamic modes of finite wavelength. Also, there is no contribution from the homo- 
geneous hydrodynamics (that for the invariants) due to the orthogonal projection (1 — P). 
The appearance of this projection is an essential self-consistency of the analysis, and occurs 
as well for normal fluids. The expression fl6.16p is only formal and the actual limit should be 
taken in the weak sense only after (16. 8p has been used to define average properties. A tech- 
nical complication is the occurrence of periodic time dependence, the Poincare recurrence 
time. This can be removed by considering the thermodynamic limit of V — > cxo, ^ oo at 
constant N/V . Therefore, averages using the normal solution to the Liouville equation are 
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understood as having the thermodynamic hmit followed by the long time limit at constant 
(a«(r);t). 

An alternative equivalent form results from performing the integral in (16.161) using the 
explicit form (16.101) and the property fl221) of the Appendix 

p„ (r, (a„ (r) ;it)) = po (r, (^^(r); t)) 

+ hm (1 - V) (e-(^^+^")*«) (r, (a„(r); t)) ■ V (a^ (r) ; t) . 

(6.17) 

The decay time for the integrand of (16.161) now becomes the time after which (16.171) reaches 
its normal form. 



A. Constitutive equations 

The exact macroscopic balance equations are given by (15. 2p . and the necessary constitu- 
tive equations are given by (15.51) and (15.61) as averages over the normal solution. These can 
be made more explicit now using the small gradient result ( 16.16^ . Since the latter is a local 
function of the fields, the constitutive equations also will be local. Furthermore, since all 
components of the gradients in (I6.17P depend on the common value (g (r) ;/!:), this can be 
eliminated through a Galilean transformation so that all properties refer to a fluid element 
at rest. Of course, the gradients of (g (r) ; t) in that fluid element are nonzero. 

Consider first the energy loss function uj 

uj{{aa (r) ;t)) = j dVpn {T, (a^ (r) ; t)) w{v) = j dVpn (r, (a^ (r) ; t)) w. (6.18) 

The coefficients of the gradient in the normal solution have translational invariance and the 
average is independent of r, except through its parameterization by {aa (r) The second 
equality takes this into account by replacing w{r) by its average w 

w = V-^ J drw{r). (6.19) 

Since u!{{aa (r) ;t)) is a scalar, fluid symmetry restricts the contributions to flrst order in 
the gradients to 

uj{{a^ (r) ; t)) = uJo{{a^ (r) ; t)) + cui((a, (r) ; t))V ■ U (r, t) . (6.20) 
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Here, the flow velocity U (r, t) of (14.41) has been used in place of the momentum density. 
The first term is the contribution from the HCS distribution 



cuo((a«(r);t)) = J dTpo{T,{aa{r);t))w. (6.21) 
The coefficient of V ■ U (r, t) is 

ui{{a^ (r) ; t)) = lim C^{to, (a, (r) ; t)) = C^{0, (a„ (r) ; t))+ lim [ dt'C^{t', (a, (r) ; t))dt' 

(6.22) 

with the correlation function defined by 

C^{t) = j dTwil~V)e-^'Mu, (6.23) 

The coefficient c^o defines an "equation of state" for the granular hydrodynamics, and 
gives the first non-trivial result of this analysis. It is similar to the pressure (given below) 
and is an inherent property of the local state of each cell, independent of the gradients 
between cells. In contrast, ui is a true transport coefficient characterizing communication 
between cells. The first equality of (16.221) provides the Helfand form for this coefficient, 
while the second equality gives the equivalent Green-Kubo form. Each has its practical 
utility, depending on the method used for its approximate evaluation. Both forms have 
proven useful for normal fluids, and further discussion is provided below. Both ujq and the 
transport coefficient ui vanish for normal fluids since they characterize collisional energy 
loss. 

The fluxes (3a of (15.61) can be determined in a similar way. As indicated in (14. 5p . only 
the rest frame flux (3'^ is required. Furthermore, since all components of the gradients 
in (I6.17P depend on the common value (g (r) ; t), this can be eliminated through a Galilean 
transformation so that all properties refer to a fluid element at rest. Of course, the gradients 
of (g (r) ; t) in that fluid element are nonzero. The component l3[ is the rest frame mass flux 
which is expected to vanish in order to give the continuity equation. This follows from the 
fact that h[ (r) is the microscopic momentum density 

P[= /cirg(r)po(r,(a4r);t))+ lim / dFg (r) (1 - P) (•■•)= (6.25) 
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The first term vanishes since (g (r) ; t) = in the rest frame, and the second term vanishes 
since (1 — V) projects orthogonal to the mass, energy, and momentum. Thus, the expected 
continuity equation is verified. 

The fluxes and jS'^ for a = 3,4,5 are the rest frame energy and momentum fluxes. 
The energy flux transforms like a vector and therefore fluid symmetry (translational and 
rotational invariance) requires that it can depend only on gradients of scalars 

(^2 = -A ((a«(r); t)) VT(r, t) - ((a,(r); t)) V (m (r) ; t) . (6.26) 

To make the connection with Fourier's law for an atomic fluid, a temperature T(r, t) has 
been introduced through the definition 

(e(r);t) = eo((m(r);t),r(r,t)). (6.27) 

For an atomic fluid the function eQ{{m{r) ]t) ,T{r,t)) is chosen to be the thermodynamic 
internal energy density. As there is no thermodynamics for a granular fluid this func- 
tion is arbitrary and simply constitutes a change of variables from (e (r) ; t) , (m (r) ; t) to 
T(r,t)), (m (r) ; t) . In this form fl6.26p is a generalization of Fourier's law where A is the 



thermal conductivity 



271] . The contribution from the gradient of the mass density is new to 



granular fluids [fi = for atomic fluids). These coefficients are given by 
A ((a„(r); t)) = lim Cx{to, (a« (r) ; t)) = Ca(0, (a, (r) ; t)) + lim [° dt'Cx{t', (a, (r) ; t))dt' 



to — 5-00 to — ^oo 







(6.28) 

fi ((a,(r); t)) = lim C^(to, (a, (r) ; t)) = C^{0, (a, (r) ; t)) + lim /" ° dt,C^{t', (a, (r) ; t))dt' 

to^co to— +00 

(6.29) 



with the correlation functions 



Cx{t) = dr(3'2 ■ (1 - P) e-(^+^-)*MT, 



(6.30) 



+ e-r^^^-> ( -^)^ iMr).) M. ) , (6.31) 



d (m (r) ; t) K22 J Ocq 
5T (r', t) 



Mt = I dr'r' ( J^""' ^ ) . (6.32) 



S{aa\t)=Q 
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. / dr'r' ( ^f°; ) . (6.33) 

Finally, the set of vectors P'^ for a = 3, 4, 5 define the pressure tensor P'^ Pij. Fluid sym- 
metry then determines that it can couple only to the momentum gradients, or equivalently 
the flow velocity gradients, in the form 

Pij = P ((aa(r); t)) % - r] ((a„(r); t)) (^d^Uj (r, t) + d,Ui (r, t) - ^5,, V ■ U (r, t)^ 

-/€((a„(r);t))5i,V-U(r,t). (6.34) 

The scalar function p ((aQ,(r); t)) is the pressure, now identified as 

p ((a„(r); t)) = [ dTpo (F, (a„ (r) ; t)) (6.35) 



The transport coefficients in fl6.34p are the shear viscosity rj {{aa{r); t)) and the bulk viscosity 
K((a„(r);t)) given by 

7] ((a,(r); t)) = lim C^(to, (r) ; t)) = C^(0, (a„ (r) ; t)) + lim [ ° dt>Cr,{t', (a, (r) ; t))dt' 

to^co to— ♦CO Jq 

(6.36) 

((a,(r); t)) = lim ^(to, (a« (r) ; t)) = ^(0, (a, (r) ; t)) + lim /" " dt^C^it', (a, (r) ; t))rft' 

(6.37) 

with the correlation functions 

Cr,{t) = j dVP'^y ■ (1 - P) e-^*M„ (6.38) 

a(t) = y c?F/3^, ■ (1 - P) e-^*M«, (6.39) 

This completes the formal derivation of the constitutive equations leading to the nonlinear 
Navier-Stokes equations, including expressions for the cooling rate, energy flux, and pressure 
tensor including contributions up through first order in the gradients of the hydrodynamic 
fields. These expressions are functions of the hydrodynamic fields to be determined by their 
detailed many-body analysis of the correlation functions. 
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B. Green-Kubo expressions 



To contrast the results here with those for an atomic fluid, it is instructive to focus on the 
Green-Kubo forms for the transport coefficients [s]. These are given by the second equalities 
of 06.221), (16.281) . 06.291) . 06.36p . and 06.37P : the first equalities are the corresponding Helfand 



forms 



28|. For an atomic fiuids there is no counter part to uJi and /i. However, there 



are Green-Kubo expressions for the thermal conductivity and the two viscosities. For the 
discussion here only the thermal conductivity is considered, whose Green-Kubo expression 
is 

A ((a,(r); t)) = Ca(0, (a„ (r) ; t)) + lim / ° dt'C^{t\ (a„ (r) ; t))dt' (6.42) 
dtCx{t) = ^-j dT!3'^ ■ (1 - P) e-(^+^-)*TA, (6.43) 

Te = -(Z+ir22)Me (6.44) 

In contrast the thermal conductivity for an atomic fiuid is 

A ((a„(r); t)) ^ lim / " dt^Cxit', (a, (r) ; t))dt' (6.45) 

dtCx{t) - ^ 1 ■ (1 - P) e-^*/5^pe, (6.46) 

In this last expression pe is the equilibrium Gibbs ensemble, and it is understood that /Jg is 
the microscopic expression for the energy fiux for a dynamics with conservative forces, and 
L is the Liouville generator for the corresponding dynamics. 

There are several similarities and differences between the granular and atomic fluid ex- 
pressions jl^. The latter is the time integral of a energy flux - energy flux equilibrium time 
correlation function. The granular fluid is similar, with one of the fluxes the same but the 
other flux is generated from the local HCS state. Also, the generator for the dynamics in the 
granular case has two additional effects, L replaced by £ + ^^225 to represents homogeneous 
cooling of the reference state and its homogeneous response to perturbations. The projection 
orthogonal to the invariants of each dynamics {1 —V) occurs in both cases as a necessary 
condition for the long time limit of the time integral, and the corresponding existence of the 
normal state. Finally, the contribution from Ca(0, {a^ (r) ;t)) vanishes for a normal fluid 
but is non-zero for the granular fluid due to the non-conservative forces. 
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C. Navier-Stokes hydrodynamic equations 



In closing this section it is appropriate to record the results of substituting the Navier- 
Stokes constitutive equations, valid to first order in the gradients, into the exact macroscopic 
balance equations. This defines the Navier-Stokes hydrodynamic equations for a granular 
fluid 

Dtm + mVr ■ U = (6.47) 
Aeo + cuo + (^p + cji + Qr] - V ■ V ■ U 
- V {daUfs + df^Ua) d^Up - V ■ (AVT + /iVm) = 0, (6.48) 



(6.49) 



For simplicity of notation, m = {m (r) ; t) in these equations. They are a set of five nonlinear 
partial differential equations for the variables m, eo, and U. They are a closed set of equations 
since uq, p, and the transport coefficients Ui, A, /i, r], and k a defined as functions of these 
variables. These definitions for the constitutive equations are the primary accomplishment 
of the statistical mechanical basis for the hydrodynamic equations. The form of fl6.36l) - (l6.37p 
could have been guessed from the outset based on the macroscopic balance equations and 
fluid symmetry. The underlying basis in the microdynamics of the particles provides the 
necessary details for how the parameters of these equations must depend on the fields. Here 
only the formal definitions have been identified. It is only the first half of the problem of 
completing these equations, as the evaluation of these definitions poses a difficult many body 
problem. Still, without this first half, the starting point for that detailed analysis would not 
be possible. This the case for atomic fluids as well. 



VII. FUTURE DIRECTIONS 



The objective here has been to formulate the basis for a macroscopic description of gran- 
ular fluids using the fundamental principles of nonequilibrium statistical mechanics. The 
analysis presented follows that for an atomic fluid. First, the exact macroscopic balance 
equations are identified. Next, their closure is linked to the concept of a normal state and 
corresponding normal solution to the Liouville equation. This defines the domain of hydro- 
dynamics in its most general sense, both for atomic and granular fluids. The construction 
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of a normal solution is quite difficult in general, but can be accomplished for states with 
small gradients relative to locally homogeneous conditions. This gives the Navier-Stokes 
approximation described here. 

Navier-Stokes hydrodynamics is applicable for most common states of atomic fluids, while 
deviations occur primarily for more complex polymeric molecular fluids. The latter have 
rheological properties corresponding to larger gradients relative to additional microscopic 
length and time scales. The construction of normal states in these cases is more difficult 
and is still at the semi-phenomenological stage [29i| . Granular fluids provide a new motivation 
for renewed efforts to describe these more complex normal states. The reason is that even 
structurally simple granular fluids composed of spherically symmetric particles can exhibit 
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3lj. This 



rheology and other phenomena beyond the Navier-Stokes domain of validity 
is due to the cooling rate in the energy balance equation which provides a new internal 
time scale, that can set the size of hydrodynamic gradients beyond any control through 
boundary conditions. For example, new steady states are possible for granular fluids due to 
the balance of this internal cooling with external forcing. In many cases this implies that the 
hydrodynamic description required is beyond the Navier-Stokes domain. The understanding 
of constitutive equations in these cases is poor at this point. It is hoped that the formal 
structure described here will provide the appropriate basis for studies of these problems. 

The context of hydrodynamics depends on the formation of a normal state from more 
complex conditions. Above this has been described qualitatively as a two stage process of 
rapid velocity relaxation in each cell to a state near the local HCS, followed by hydrodynamic 
relaxation through exchange of mass, energy, and momentum between the cells on a longer 
time scale. This separation of microscopic and hydrodynamic time scales is essential to 
the dominance of the hydrodynamic excitations over all others at large space and time 
scales. It is justified for atomic fluids since the hydrodynamic times are determined by the 
wavelength of the phenomena studied. As the system approaches homogeneity, these time 
scales become much larger than the microscopic excitations and hydrodynamics prevails at 
large times. However, there is an additional hydrodynamic time scale for granular fluids, the 
cooling rate, which is not set by the wavelength alone. It would seem that this additional 
time scale must be large as well, implying a weak cooling rate. This condition is too strong. 
What matters is the rate of the approach to the homogeneous state, not any dynamics of 
that flnal state. In the above derivation of hydrodynamics the flnal form for the solution to 
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the Liouville equation, eq. fl6.16p or (16 .17^ . has a dynamics generated by IC + rather 
than simply that for the trajectories L. This is significant since the former has the additional 
compensation for the cooling and for the homogeneous perturbations of that cooling. Hence 
the approach to the time dependent normal state is determined only by the remaining non- 
hydrodynamic relaxation. The time scale for relaxation to the normal state is independent 
of the hydrodynamic time scales of that normal state. Quantitative verification of these 
concepts is another important future direction for research on a hydrodynamic description 
for granular fiuids. 
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Appendix - Gradient expansion 

In this Appendix the Liouville equation in the form (16. 6p is written to first order in the 
gradients and solved. Also the invariants of the associated dynamics are identified. 
Consider first the right side of (16.61) which can be written equivalently as 



dr' . ^. {V ■ (b„ (r) ; t) + 5^, {w (r) ; t)} - Lp, 



I 

= - [ dr' {La^ (r') ; t) - Lpoe 

J 5 {aa (r') ; t) 

The first equality follows from (14.21) and (15. ip . The first two terms are determined by the 
local HCS which can be expanded to first order in the gradients 

Poe = Po {{aa (r) ; t)) + / dv' ( ) ((a„ (r) ; t) - (a„ (r) ; t)) + ■ ■ 

= Po ((a, (r) ; t)) + (r, (a, (r) ; t)) ■ V {ap (r) ; t) + •• (2) 

The functional derivatives are 

^Poe \ . f dpo{{a^{r);t)) (r, (a„ (r) ; t)) 



'^(«a(r');i)/5(a,;t>=o V 5(a„(r);t) 9 (a« (r) ; t) 

+m,(r,(a,(r);t)) • V5(r'-r) + ■■ (3) 



Here 



m, (r, (a„ (r) ; t)) ^ / dr' (jr^f^) (r'-r) , (4) 

and Po ((o-a (i") is the actual HCS with its global density, energy, and momentum evalu- 
ated at the common values {aa (r) It follows from (16.41) that the averages of Oq, (r) for 
p,Poe, and po are all the same. This in turn gives 

dVaa (r) ni/3 = = / dVaa (r) A. (5) 
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With these results and the fact that A is of first order in the gradients, ([T]) to first order 
in the gradients becomes 

f ^^ ^/ A ■ ' ^) + '^"2 {w (r) ■,t)}- Lpoe 

J b {tta (r'j ; t) 

^ Cpo -{l-V) [IC + K^)^^ uif, ■ V (a,; t) + VCA (6) 
The matrix is the transpose of K 

duj{{aa (r);t)) 
- ^'^^ d{a,{r);t) ' 

and J is the unit matrix. The generator C is the same as that of fl3.12p with u 
i^o((oa (r) ; t)) for the HCS evaluated at the common values (r) ; t) 

C = -UJo{{aa{T);t))d(^e{r);t)+ L. (8) 

Finally, V is the projection operator 

= ^/ / rfra„ (r) X. (9) 



9 (a„ (r) ; t) 

The first term of ([6]) vanishes by definition of the HCS, po, confirming that the right side of 
the Liouville equation fl6.6p is of first order in the gradients. 
At this point, the Liouville equation (16.61) becomes 

dtA - [ dr'— 4^-— a;o(K (r) ; ^)) + VLA = {l-V)K-V (a„; t) , (10) 

T^^-{IC + K^)^^m, (11) 

This equation is still exact up through contributions of first order in the gradients. It has 
solutions of the form 

A (r, t I (a, (r) ; t)) = (F, t, (a„ (r) ; t)) ■ V (a, (r) ; t) , (12) 

Substitution into (|TOl) gives the corresponding equation for 

dtG. + (1 - P) [IC + K^)^^ Gp = {l- V) n, (13) 

with the solution 

(r, t, (a„ (r) ; t)) = f dt' (e-(i-^)(^^+^n*') (1 _ V) T,. (14) 



28 



It is possible to add to (fTOj) an arbitrary solution to the homogeneous equation corresponding 
to f|T3|) . As described in the text, this represents the dynamics of the first stage of rapid 
velocity relaxation to the local HCS. The interest here is in the second stage where possible 
formation of a normal solution occurs. Hence, it is simpler to choose that stage for initial 
conditions (initial local HCS). 

Define the derivatives of the HCS by 



Then differentiate the equation for po 

d 



-Cpo = 0, (16) 



d {a/3 (r) ; t) 
to get 

{ICt + K'')^^^, = 0. (17) 

Since [ICt + K^) is the generator for the dynamics in (fT^ this shows that \l//3 are the 
invariants of that dynamics. 

The projection operator V in (fT8|) acts only on phase functions with translational invari- 
ance. In that case iQ simplifies to 



VX = ^p j dVApX, = V-^ j draa (r) . (18) 

The first equality of ([5]) becomes Pm^j = 0. This in turn gives 

m,= (l-P)m,= (l-P)M,, = / dr' f — %— ^) r'. (19) 

Then (1 —V)T'^ simplifies to 

{l-V)T'^ = -{l- V) {IC + K^)^^ (1 -V)M^^{1- V) T„ (20) 

with 

T^ = -{IZ+K^)^^m,. (21) 
Use has been made of the identity 

{1-V){IC + K^)V = Q. (22) 

This same identity leads to a simplification of the dynamics in (JT^ 

^_(i_P)(,S+,,T),, _ ^ _ g-(/£+/<-)i' (1 _ p) . (23) 
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In summary, the solution to the Liouville equation to first order in the gradients is 

p{r,t\{aa;t)) = po(r,(a„(r);i)) + (l-P)(M^(r,(a„(r);i)) 

+ 1^ dt' (e-(^^+^")*') (1 - V) Tfs (r, {a^(r);t))^ ■ V (a. (r) ; t) 

(24) 
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